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ESTIMATION OF THE TIME OF CHANGE IN PANEL DATA 


LAJOS HORVATH, MARIE HUSKOVA, GREGORY RICE, AND JIA WANG 


Abstract. We consider the problem of estimating the common time of a change in the 
mean parameters of panel data when dependence is allowed between the panels in the form 
of a common factor. A CUSUM type estimator is proposed, and we establish first and 
second order asymptotics that can be used to derive consistent confidence intervals for the 
time of change. Our results improve upon existing theory in two primary directions. Firstly, 
the conditions we impose on the model errors only pertain to the order of their long run 
moments, and hence our results hold for nearly all stationary time series models of interest, 
including nonlinear time series like the ARCH and GARCH processes. Secondly, we study 
how the asymptotic distribution and norming sequences of the estimator depend on the 
magnitude of the changes in each panel and the common factor loadings. The performance 
of our results in finite samples is demonstrated with a Monte Carlo simulation study, and we 
consider applications to two real data sets: the exchange rates of 23 currencies with respect 
to the US dollar, and the GDP per capita in 113 countries. 


1. INTRODUCTION 

In this paper, we consider the problem of estimating the time of a change in the mean 
present in panel data in which their are N panels comprised of time series data of length T. 
A common structural break in panel data is a quite natural occurrence. For example, if each 
panel represents the exchange rate of a currency with respect to US dollars, then a crisis in 
the US would be expected to simultaneously affect each panel. Similar phenomena may be 
produced my governmental policy changes, the introduction of a new technology, etc., and 
in these cases it is of interest to estimate the time at which such occurrences are manifested 
in sample data. The theory of change point analysis has been extensively developed to study 
problems of this nature, see Csorgo and Horvath (1997), Brodsky and Darkhovskii (2002), 
and Aue and Horvath (2012) for reviews of the held. 

Classical methods in change point analysis consider univariate and multivariate data of 
a hxed dimension. In many panel data examples, however, the number of panels N is 
comparable in size to the length of the series T. In these cases asymptotics as T remains hxed 
and N tends to inhnity, or as N and T jointly tend to inhnity, are more appropriate. Although 
in principle one could detect the common change present in each panel by examining a single 
series, an analysis that utilizes all available panels should provide improved detection and 
estimation. 

The literature on structural breaks in panel data has grown considerably in the last two 
decades. We refer to Arellano (2003), Hsiao (2003) and Baltagi (2013) for surveys of several 
panel data models and their applications to econometrics and hnance. The early founda¬ 
tions for estimating structural changes in panel data were developed in Joseph and Wolfson 
(1992,1993), and many aspects of the problem have now seen at least some consideration; 
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Li et al. (2014) and Qian and Su (2014) consider multiple structural breaks in panel data, 
and Kao et al. (2014) considers break testing under cointegration. 

The works of Bai (2010), Kim (2011, 2014), and Horvath and Huskova (2012) are the most 
closely related to the present paper. Bai (2010) considers the problem of estimating a 
common break in the means of panel data that do not exhibit cross sectional dependence. A 
least squares estimator is proposed that is shown to be consistent when N tends to inhnity, 
and its asymptotic properties are derived as N and T jointly tend to inhnity. Kim (2011) 
considers the least squares estimator of Bai (2010) with cross sectional dependence modeled 
by a common factor, and expands the test to detect a change in the slope of a linear trend in 
the mean component. Horvath and Huskova (2012) study testing for the presence of a change 
using a CUSUM estimator, and also assuming the presence of a common factor. In each of 
these papers, asymptotics are derived assuming the model errors are linear processes, and 
that the rates of divergence relative to N and T of the size of the changes and the magnitudes 
of the factor loadings are hxed. 

In this paper, we expand the existing theory in two primary directions. We derive second 
order asymptotics for the CUSUM change point estimator assuming only an order condition 
on the long run moments. This extends the asymptotic theory of change point estimation 
to a wide variety of error processes, notably many nonlinear time series examples like the 
ARCH and GARCH processes. We also show explicitly how the asymptotic distribution and 
norming sequences of the estimator depend on the magnitude of the changes in each panel and 
the common factor loadings. This allows for the computation of the limit distribution under 
several conceivable rates of divergence for the magnitudes of changes and factor loadings. 
The remainder of the paper is organized as follows. In Section we present our assumptions 
and the main results of the paper. Section [^contains examples of error processes that satisfy 
the assumptions of Section]^ Estimators for the norming sequences that appear in the results 
of Section are developed and studied in Section The implementation of the results of 
the paper as well as a Monte Carlo simulation study and data applications are detailed in 
Section The proofs of all results are contained in Appendix 

2. Assumptions, and main results 
We consider the panel data model 

Xi^t = 1^1 +5il{t > to} +'yir]t +ei^t, ^<i<N,l<t<T, (2.1) 

where the idiosyncratic errors e'^s have mean zero, r;* denotes the common factor with 
loadings y*, 1 < i < N, and hj denotes the change in the mean of panel i that occurs at the 
common, and unknown, change point to- 

Assumption 2.1. 

(i) The sequences —oo < t < cx)}, 1 < i < N are independent, and 


(ii) {rjt, —oo < t < cxd} and —oo < t < cx)} 1 <i < N are stationary. 

According to Assumption 2.1[ i), the only source of dependence between the panels is the 
common factor 7]t- The idiosyncratic errors form a stationary time series, similarly to the 
assumption in Bai (2010) and Kao et al. (2012, 2013). Throughout this paper 5i and y,, 
1 < z < A, are allowed to depend on N and T. For the sake of simplicity, we consider the 
case when y^ G R, but our results could be extended to the more general case of a vector 
valued factor loading and common factor. 
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Assumption 2.2. The time of change in the mean to satishes 

to = [T6\ with some 0 < 6^ < 1. 


Assumption |2.2 is standard in change point analysis, and corresponds with the assumptions 


of Bai (2010), Kim (2011,2014), and Horvath and Huskova (2012). It is of interest in some 
econometric applications to allow for 9 to depend on N and T and tend to the end points 
0 or 1 at a certain rate; see Andrews (2003) and, in the panel data setting, Qian and Su 
(2014). The consideration of this problem in generality for our estimator is not a goal of the 
present paper, and requires a thorough study. 

Our estimator for to is dehned as the location of the maximum of the sum of the CUSUM 
processes across the panels: 



i=l 


where 




S = 1 


The estimator of Bai (2010) is 



( 2 . 2 ) 


2=1 


which is the maximum likelihood estimator for to assuming that the panels are independent 
and normally distributed with the same variance, while tN,T maximizes the weighted log 
likelihood. 

We impose only conditions on the long run moments of the error processes for our asymptotic 
results. The long run moments of the errors in panel i are dehned by 


= E , I <i <N, 


S=1 

and we assume that they satisfy the following conditions: 


Assumption 2.3. 

(i) There exists Ci,! < i < N such that 


m^ sup -f/i, 2 (t) - cr^ 

i<i<N t 


where Ci < ai < C 2 for all 1 < i < N with some 0 < Ci < (72 < 00 , and 

(ii) 



Additionally, we must assume an analogous condition on the common factors: 


Assumption 2.4. With some k > 2. 


Erjt = 0, E = t + o{t), and E rjg 



S=1 / S = 1 


ast^oo. 
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Assumptions 2.3 and 2.4 do not assume any specific structure on the error terms, in contrast 
to the structural break literature with panel data to date. We provide several examples 
in Section including linear and nonlinear time series, martingales and mixing sequences, 
where Assumptions 2.3 and 2.4 are satished. 


The size of the changes and the correlation between the panels will play a crucial role in the 
asymptotic distribution of the estimator, and these quantities will be measured by 

N N N 

^N,T = ^ ^ Sjv,T = ^ Si'Ji- 

2=1 2=1 2=1 

The limit results below are proven when min(A^, T) —)■ oo. 

Assumption 2.5. As min(A, T) —)■ oo, 


TA 


lb 

and 


N,T 


N 


—)■ oo, 


11 


N,T 




0 . 


Assumption 2.5 means that the sizes of all changes cannot be too small and that the factor 


loadings cannot be much larger than the sample size and the size of the changes. Bai (2010) 
assumes that Ajv,t/A^ converges to a positive limit while under the assumptions of Kim 
(2011), the common factor dominates. A primary goal of our paper is to show how the 
relationship between the loadings and the sizes of the changes affect the limit distribution 
of the time of change estimator. 

Our hrst result pertains to the asymptotic distribution of tN,T when An^t is large. 


Theorem 2.1. If Assumptions 2.1- 2.5 


and 


as N,T —)■ oo, then we have that 


At oo 

(2.3) 

f'" = 

^N,T 

(2.4) 

P = ^o} —t 1. 

(2.5) 


The assumption in (2.4) may seem somewhat restrictive since it rules out the example of 


hxed break sizes and factor loadings. We note that due to the result on page 635 of Horvath 
and Huskova (2012), when the factor loadings are hxed the CUSUM test for the presence of 
a change point will reject with probability tending to one regardless of if a change exists or 


not, and so something along the lines of (2.4) must be assumed for the CUSUM estimator 
of the time of change to be consistent. 


Remark 2.1. Assume that T is hxed. If Assumptions 2.1 -2.4 are satished and T ^n,t 


0, and An,t/N —)■ oo, then (2.5) holds. 
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Remark 2.2. In order to establish the consistency of Bai’s (2010) estimator in (2.2) for hxed 
T and under our assumptions we must assume in addition that for each and r]t are 

uncorrelated random variables, and that {e^, 0 < t < oo , I < i < N} and {rjt, 0 < t < oo} 
are independent. If in addition to Assumptions 2.1-2.4, EtJq < oo, —>■ oo and 

r isf^T/^N,T —t 0 hold, then we have that 

lim P{t%j, = to} = 1. (2.6) 

N^oo ’ 


We provide a proof of Remark 2^ in Appendix 


The main difference between Remarks 
cx). Remark 


2.2 


2.1 


and 2.2 is in the assumptions Ajv,t/A —)■ oo and 


allows smaller changes to establish consistency but much 
stronger assumptions on the sequences ei^t, I < i < N and rjt. If we cannot assume that 
Gi,t, 1 < i < N and rjt are sequences of uncorrelated random variables and the independence 


of {vt^t > 0} and {ci^tP > 0,1 < i < N}, then (2.6) can be proven under conditions of 
Remark 2.1 In this case (2.5) and (2.6) can hold only if —>■ oo when T is hxed. 


We now turn to the asymptotic distribution of AAr,T(^A,T — ^o) when (2.3) does not hold, i.e. 
the sizes of the changes are small or occur in only a few panels: 

Assumption 2.6. 

(i) Atv.t = 0(1), 


(ii) TA7 v,t( 1 + log(T/AAr^T)) —)■ oo, where k is dehned in Assumption 2.4, and 


(iii) “^PAn^t/N^P —)■ oo, where k is dehned in Assumption 2.3 (ii) 


2.5 


we have that TAj^^t —t oo, so Assumption 


2.6 


By Assumption 

oo. If N —>■ oo, i.e. the number of the panels is large, then Assumption 


ii) holds if, A^T(logT) 
follows 


from Assumption 2.5 However, Assumption |2.5[ii) also holds if the number of panels is 


2.6 


111 


relatively small with respect to the length of the panels and the sizes of the changes. 


Next we introduce an assumption that is a companion to Assumption 2.3 

Assumption 2.7. 


N 


(i) — \Eeiflei^t\ = 0{t ’’) with some r > 2, and. 


2=1 


(ii) with some r > 2. 


Assumption 2.7 requires an upper bound for the average correlation of the errors in the 
panels and a uniformity condition that augments Assumption 2.3[ii). 


Our hrst result in this direction covers the case when the sizes of the changes are small and 
the effect of the correlation between the panels is negligible or moderate. We measure the 
dependence between the panels with respect to the sizes of the changes by 


5 = 


min 


lim — 

{N,T)^oo 


N,T 

lj2- 

N,T 
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To describe the limit distribution of t 7 v,T we need to introduce a drift function 

f (1 — 6')|m|, if m < 0 

ge{u) = 

[ 9u, if M > 0 

and an asymptotic variance term 

= lim 


N 


N,T^oo 




V?, 


2 = 1 


We note that by Assumption |2.3Ki) we get that Ci < a < C 2 - Let 


/ t 




S=1 


if t = 1,2,3,. 

Ui{t) = •^ 0, if f = 0 

-1 


if ^ f 5 2, 3,.. 


S = t 


The function u{s,t) is the asymptotic covariance of i-e. for all integers s and t 


N 


N 


u(s,t) = \im E j y^^Sjhlj{s) I j '^6iUi{t) | . 
i=i 


N^OO 
tTV 


. 2=1 


We note that u{s,t) = hm 7 v-).oo X]i=i EUi{s)Ui{t). It follows from Assumpti 


ion 


2.3 


(2.7) 


i) that 


the covariance function u(s, t) is hnite for all integers s and t. This function only appears in 
Theorem 2.2 below when Ajv,r is above some positive bound for all N and T. In this case 
u(t, t) > 0, if t 7^ 0. 

The next result considers the case when the common factors are negligible. 


Theorem 2.2. We assume that Assumptions 2.1-2.1 hold, 


N 


^N,T l'^*r 0, 


2=1 


where r is defined in Assumption \2.1^ ii), and 

5 = 0. 

(a) If 

Ai.N,T —t 0, 

then we have 


Il^N,T(tN,T — to) V X , X. 

- argmax„ {W[u)- geW)] , 




where W{u), —oo < u < oo is a two-sided Wiener process. 

(b) If 

I^N,T —t h G (0, oo), 

then we have 

tN,T -to -> argmax,{0(«) - 5g»(«)}, 

where t — 0, ±1, ±2, ... is Gaussian with E0(t) — 0 and E0(s)0(t) — 11 ( 5 , t). 


( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 
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Remark 2.3. Since the proofs of Theorem 2.2 and the results to follow depend on nor 


mal approximations for the sums fhe independence of 

> 0} on z could be relaxed as it is pointed out by Bai (2010) and Kim (2011). This 
would be an important consideration, for example, if the panels are indexed by locations, i.e. 
z = i describes the location of the panel. In this case a spatial structure could be assumed 
on the errors. If Assumption |2.lK i) is replaced by a weak dependence or spatial assumption, 
the norming constants would change in our limit theorems. 


N 


NT 0 with some r > 2 does not hold, then the limit in (2.13) 


Remark 2.4. If A 

might not be the argmax of a normal process with a drift 


For example, if (5i = 1 and 5^ = 0 


for all z > 2, then the limit in (2.13) is determined by the error terms ei^t, —oo <t < oo. 


Remark 2.5. The distribution of argmax„{hF(M) — geiu)} is known explicitly. Its density 
was derived by Ferger (1994) from Corollary 4 of Bhattacharya and Brockwell (1976) (cf. 
Csorgd and Horvath (1997, p. 177)). 


Remark 2.6. If (2.12) holds, 7 * = 0 for all 1 < z < A and EeiflCi^t = 0 for all 1 < z < A 
and t 7 ^ 0, then we get the analogue of Theorem 4.2 of Bai (2010). In this case 0(s) is a 


Wiener process on integers, so the main difference between the limits in (2.11) and (2.13) is 
that the argmax is computed on the real line or on integers. 

So far the common factor was part of the error term and it had a negligible contribution to 
the limit distribution in Theorem 2.2 However, it has been observed in testing for changes 


in panel data that the effect of strong correlation between panels might make standard 
statistical procedures invalid (cf. Horvath and Huskova (2012)). Bai and Zhang (2010) 
provide some economic examples when the common factor between the panels is strong so 
its effect must be built into testing and estimating procedures. Our next result covers the 
case when the order of the correlation between the panels and the sizes of the changes are 
essentially the same. Since the contribution of the to the limit will not be negligible we 
need to specify the relation between the errors and the common factors: 

Assumption 2.8. The sequences {rjt, —oo < t < oo} and —C )0 <t<oo,l<z< 
A} are independent. 


Similarly to Ui(t), we introduce 






S = 1 


V(t) = 0, 


if t = 1,2,3,. 
if f = 0 


if f =-1,-2,-3,..., 


S=t 


which will be part of the limit distribution when (2.12) holds or is proportional to 

An,t- Ih all the other cases we assume the asymptotic normality of V{t): 


Assumption 2.9. 


[Tu\ 


Y, m 


D[0,1] 


IK(zz), where IF is a Wiener process. 


z=i 
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Theorem 2.3. We assume that Assumptions 2.1-2.9, and (2.8) hold, and 


(a) If (2.10) holds, then we have 

A^N,TitN,T — to) 


0 < |s| < cx). 


V 


-)■ argmax„ {W{u)- ge{u)} , 


+ 5 ^ 

where W{u), —oo < u < oc is a two-sided Wiener process, 
(b) If (2.12) holds, then we have 


(2.14) 


(2.15) 


iN,T -to 4 argmaxj{(S(t) + sd^^‘^V{t) - Dgeit)}, 


(2.16) 


where <S(t),t = 0, ±1, ±2,... is the Gaussian process defined in Theorem 2.2, independent 
ofV{t),t = 0,±l,±2,... 


The effect of correlation between the panels is demonstrated in Theorem 2.3 The limit 


distribution in (2.15) remained the same as in (2.11) but the variance of the estimator 
increased by s^. The effect of the common factor is more transparent in (2.16) since an 


additional term appears in the limit which depends on the distribution of the common 
factors. 


Theorem 2.4. If Assumptions 2.T2.9, and (2.8) hold, and 


then we have 


| 5 | = OO, 


A2 

v^(tA,r-to) argmax„{lT(M)-5(0 (m)}, 

^N,T 


(2.17) 

(2.18) 


where W{u), —oo < u < oo is a two-sided Wiener process. 


Theorem 12.41 covers the case when the limit distribution of the estimator for the time of 


change is completely determined by the common factors. The limit distribution in (2.18) is 


the same as in (2.11) and (2.15) but the rate of convergence is much slower. The effect of 


having several panels with changes is overshadowed by the strong influence of the common 
factor. For further results when the common factors dominate the limit distribution we refer 
to Kim (2011). 


3. Examples 


In this section, we study some examples of error processes that satisfy the assumptions of 
Section We restrict our attention to establishing Assumption 2^ for examples of possible 
model error sequences Cj^ds, but the same sequences could be used for the common factors 
as well. 


Example 3.1. Let —oo < t < oo be independent, identically distributed random vari¬ 
ables with Ecifi = 0, Ee^Q = af and Elej^oT < Due to independence we have that 

Ui, 2 {t) = taj for alH = 0,1, 2 ... 


(3.1) 
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By the Rosenthal inequality (cf. Petrov (1995, p. 59)) we obtain for all t > 1 that 

U,,.{t) < C {tS|e,,or + < Ce'^{E\ei^o\^ + af}, 

where C is an absolute constant, depending on k > 2 only. If the error terms in each panel 
are independent and identically distributed, then, assuming Ci < cxj < C 2 for all 1 < z < A^, 
Assumption 2.3Ki) holds; Assumption 2.3[ ii) is satished if 

N 


N 


ol'")^ = 0(1) with some a > 4. 


2=1 


If maxi<i<Ar h^|ei,or — ^0 with some r > 2 and Oq > 0, then Assumption 
fulhlled. 


(3.2) 

2^ii) is also 


ARMA processes are very often used in classical time series analysis and our next example 
shows that stationary ARMA processes satisfy the basic assumptions of the hrst section. 
We consider the more general case of linear processes, which are investigated by Bai (2010), 
Kim (2011), and Horvath and Huskova (2012). 


Example 3.2. We assume that —00 < t < cx)} are independent and identically dis¬ 
tributed random variables with Eei^ = 0 and E\eifl\^ < 00 with some k > 4. The error 
terms form a linear process given by 




00 

£=0 


where sup^^ ^ < Q with some Q > 0 and Oj > 0. By the Phillips and Solo (1992) 

representation we get 

t t 00 / t \ 

s=l s=l j=—oo \fc=l / 


where Q = ^ 0 ’*^*-0 = Ci,o - C'i,Ciy = Ci^i, if £ > 1 and = 0, if £ < -1. 


Minkowski’s inequality anc 
CfEef^ in Assumption 
as t —)■ 00 . 


2.3 


the discussion in Example 3.1 yield that we need to choose af = 
i) and we also have Ui^nit) = and |Eei,oej,t| = 0(t“^““‘), 


The next example assumes a martingale structure of the errors. 

Example 3.3. Let us assume that —00 < t < cxd} is a stationary orthogonal martingale 
difference sequence with respect to some hltration with Eci^ = 0, Eej^ = af and Elej^oT < 
CX). Then (3.1) as well as Assumption 2.7 hold. By Li (2003) we also have 

Ui^nit) < E\eifl\^ with some constant C depending only on k, 

completing the proof of (3.2). Under assumption maxi<j<7v Elej^oT — ^0 with some r > 2 


and Co > 0, we obtain Assumption 2.7 


Since the early 1980’s, ARCH, GARCH processes and their various extensions have become 
extremely popular models in the analysis of macroeconomic and hnancial data. For a survey 
and detailed study of volatility models we refer to Francq and Zakoi'an (2010). The next 
example shows that a large class of volatility processes satishes the assumptions in Section 
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Example 3.4. We assume that —oo < t < 00 } are independent and identically dis¬ 
tributed random variables with Esifi = 0 and Ee'^q = 1. The error terms are dehned by 

(^i,t (3.3) 


where the volatility process > 0 is measurable with respect to the a-algebra generated 
by s < t — 1. Usually, hi^t is given by a recursion involving s < t — 1. Francq 


and Zakoian (2010) provides conditions for the existence of a stationary solution of (3.3) 
in several models and establishes their basic properties. Assuming that E\eifl\'^ < 00 with 
some A > 4, —00 < t < cx)} is a stationary orthogonal martingale satisfying the 

conditions in Example 3.3 In case of the most popular GARCH(p, q) model hi^t = u: + 


+ 'm=i > 0 , «£ > 0 ,/dj > 0,1 < £ < p, 1 < j < g. The necessary 

and sufficient condition for the existence of the higher moments in case of GARCH(1,1) is 
given in Nelson (1990). He and Terasvirta (1999), Ling and McAleer (2002) and Berkes et 
al. (2003) partially extends his results to the more general case. The existence of moments 
of augmented GARGH sequences are discussed in Garrasco and Ghen (2002) and Hermann 
(2008). 


Linear processes and the volatility models of Example |3.4| are in the class of m-decomposable 
processes. 

Example 3.5. We say the Ci^t is a Bernoulli shift if it can be written as 

= /i(G,0 G,i-l5 G,t-25 • • •) 

with some functional /*, where —oo<t<oo} are independent and identically dis¬ 
tributed random variables. The conditions of Section |2] are satished if the Bernoulli shift is 
m-decomposable, i.e. if 

00 

^(E|ej,t — < 00 with some k > 4, 

m=l 

where = /i(£j,o G,t-i, ■ • •, and the are independent 

copies of £^, 0 , independent of 1 < < A^, —00 <t < 00 . Berkes et al. (2011) proved that 
there is constant Ci such that ^ \Eeiflei^t\ < and Uj, 2 (t) = tal -|- (f), 

as f —)■ 00 , with some af. They also provide several examples for m-decomposable Bernoulli 
shifts. 


Example 3.6. There is a well developed theory of partial sums of mixing random variables 
where the long run moments play a crucial role. It has been established under various 

conditions that Ui^ 2 {t) = ierf + o{t) and = 0(G/^), as t —)■ cxd. For surveys on mixing 

processes we refer to Bradley (2007) and Dedecker et al. (2007). 


Example 3.7. We assumed in Examples |3.2| and 3.4 that the innovations Si^t, —00 < t < 00 
are independent and identically distributed. However, this assumption can be replaced with 
the less restrictive requirement that Si^t, —00 < t < 00 is a. stationary sequence. Rosenthal- 
type inequalities for sums of functionals of stationary processes are developed in Wu (2002) 
and Merlevede et al. (2006). These results can be used to establish Assumptions 2.3, and 

o 
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4. Estimation of norming sequences 


Theorems 


2.2 


2.4 


contain the limit distribution of tTvr with different normalizations to show 


the effects of the sizes of changes and the loading factors. However, in case of finite N and 
T it is impossible to check which specific condition on the growths of and Sjv,r holds, 
so we require norming sequences that would work in all possible cases. Let 

N 


^ 0 -^ 6 ^ + S 


2 

N,T- 


i=l 


Under the conditions of Theorems I2.2ll2~4l we have that 


A2 


N,T 


-Af,T 


(tN,T — to) converges in distribution. 


(4.1) 


The limit distribution in (4.1) is argmax^{lU(M) —ge{u)} except in the special cases of ( |2.13[ ) 
and (2.16). In these cases, the limit distribution is the argmax of a process defined on integers. 
The limits in (2.13) and (2.16) depend on the distributions of 1 < i < N,1 < t < T}. 
If h is close to 0 in (2.13) or (2.16) the limiting distributions can be well approximated with 
argmax^(hU {u)—g 0 {u)). Bai (2010) provides numerical evidence that argmax^jlU {u)—g 0 {u)} 
gives a reasonable approximation for the limit in (2.13) when u{s,t) = min(t,s). Hence we 
recommend that argmax„(H^(M) — g 0 {u)) can be used as the limit in (4.1|) in practice. 


The limit result in (4.1) can only be used for hypothesis testing or conhdence intervals if the 
norming factor can be consistently estimated from the sample. We estimate An^t with 

\ ^ 


2 = 1 


tN,T 






i.t 




It is more difficult to estimate Let 

N 


2=1 


UnW = ^ S,{t) - -St{T) 


rN,T{t) = -t 


T — tN^T 


and 


rN,T — rN,TitN,T) — —tN,T 

The estimator for Stv.t is defined as 

- 1 1 


T 

+ {t — t7V,T)L{t > tw.T}, 

T — 


(4.2) 


T 


(4.3) 


-A,r — 


E 

Mi<|'u|<M2 


2(M2 - Ml) ^ 4|n|f2 


N,T 


UN{tN,T + P) — UN{tN,T) 


AN,T{rN,TitN,T + t ) — 


Theorem 4.1. (i) We assume that the conditions of Theorem 2.2 or 2.3 are satisfied. If 
Ml < M 2 , Ml —>■ cx) and M 2 /T —)• 0, then 


A 


N,T 


A 


p 
—> 


N,T 


(4.4) 
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and 


■-N,T 

-7V,T 


A 1. 


(4.6) 


(ii) We assume that the conditions of Theorem 2.4\ are satisfied. If Mi < M 2 , Mi —>■ 00 and 

M 2 / vavaiT, Y,rp j^^ rp^ ^ 0 , 


then (4.4) and (4.5) hold. 


Since (4.4) and (4.5) hold, the limit result in (4.1) remains true when the norming is replaced 


with the corresponding estimators, i.e. 

T ^ 

’ -(tN,T — to) converges in distribution. 


A2 


(4.6) 


-A,T 


If the interaction between the panels is small, i.e. Y^rp/ Yl!i=i t 0, as A^, T —)■ oc, then 


we need to estimate only Ajv,r and Using df, 1 < i < N, the long run vari¬ 
ance estimators for af, a possible estimator for U ~ 

5. Simulations, and data examples 

5.1. Simulations. Using the estimators dehned in Section]^ we computed the empirical 
percentages when the variable dehned in (3.5) is below the asymptotic quantiles for several 
N and T. We considered the case when there is no interaction between the panels, i.e. 
7 i = 0 and also cases with correlations. We tried various values of to = [T^J. We observed 
that t he ap plic abili ty of the limit results presented in Section does not depend on 6. In 


Tables 
Th 


5.1 


and 

Ah 


5.2 


we used 6 = 1/2 and they are based on 1,000 repetitions. The 90 


th 


95^^^ and 99^^^ percentiles of the distribution of argmax„(lU(M) — gii 2 {u)) are 4.70, 7.69 and 
15.89, respectively. Table 5.1 illustrates that must be small to use Theorem |2.2[a) 


and the approximation imp roves when T increases. In case of larger Ajv^t we are under the 
conditions of Theorem 


2.1 


get from the asymptotics in Theorem 2.2[a). 


and the distribution of tN,T is more concentrated than we would 


Table 5.1. Empirical percentages of A^ 2 ’(^y,t — to)/Yisf^T below the asymp¬ 
totic quantiles for various values of N, T, 6i when 7 ^ = 0 for all 1 < i < A. 


N/T <5, 

90% 

95% 

99% 

25/100 

0.150 

88.7% 

94.7% 

100 % 

25/250 

0.100 

87.6% 

94.4% 

99.8% 

25/500 

0.060 

89.9% 

96.4% 

100 % 

50/100 

0.100 

89.1% 

97.1% 

100 % 

50/250 

0.070 

88 . 1 % 

95.4% 

100 % 

50/500 

0.050 

86.5% 

94.9% 

100 % 

100/100 

0.085 

86 . 1 % 

94.5% 

100 % 

100/250 

0.055 

86.9% 

94.3% 

99.9% 

100/500 

0.035 

88.5% 

96.5% 

100 % 
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Table 5.2. Empirical percentages of — to)l'^N,T below the asymp¬ 

totic quantiles for various values of iV, T, 5i when 7 * = 0.03 for all 1 < z < A. 


N/T 

90% 

95% 

99% 

25/100 

0.150 

86.9% 

94.9% 

100 % 

25/250 

0.100 

90.7% 

95.6% 

100 % 

25/500 

0.060 

89.2% 

97.2% 

100 % 

50/100 

0.100 

87.6% 

95.6% 

100 % 

50/250 

0.070 

91.0% 

96.4% 

100 % 

50/500 

0.050 

88.3% 

96.3% 

100 % 

100/100 

0.085 

89.4% 

95.9% 

100 % 

100/250 

0.055 

85.1% 

94.0% 

99.9% 

100/500 

0.035 

90.1% 

96.5% 

100 % 


Figlfre 5.1. The histogram of fioo.soo with 5i = 0.07, 7 * = 0,1 < z < 100 (left 
panel) and the histogram of ^ 50,100 with 5i = 0.1, 7 * = 0.03,1 < z < 50 (right 
panel) and the density of the limiting random variable 




Our comments also holds when interaction between the panels is allowed as illustrated by 


Table 5.2 but, due to the dependence, larger T is needed to use limit results. 


Figure [5T] shows that the density function of the limit follows the shape of the histogram of 
t]^^T closely. 

5.2. Applications. In the hrst example we consider the exchange rates between the US dol¬ 
lar and 23 other currencies. The data can be found at the website www. f ederalreserve. gov 
/releases/hlO/hist/. Figure 5.2| contains the graphs of the exchange rates between 
the United Kingdom (UK), Canada (CA), Singapore (SI), Switzerland (SW), Denmark 
(DN), Norway (NO) and Sweden (SD). In our study we used the time period 03/13/2001- 
03/11/2003 so we have N = 23 panels and each panel has T = 500 observations. Using the 
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Figure 5.2. The graphs exchange rates, 1=UK, 2=SI, 3=CA, 4=SW (left 
panel); 1=DN, 2=NO, 3=SD (right panel) with respect to the US dollar. 




Figure 5.3. The graphs of relative exchange rates, 1=UK, 2=SI, 3=CA, 
4=SW (left panel); 1=DN, 2=NO, 3=SD (right panel) with respect to the US 
dollar. 




testing method in Horvath and Huskova (2013) the no change in the mean of the panels 
null hypothesis is rejected. The estimated time of change is ^ 23,500 = 297 so the change is 



occurs in a panel with high values, this change will give a very large 6“^ compared to the 
other panels. Hence a single panel can disproportionately contribute to An,t- To overcome 
this problem we rescaled the observations in each panel with the hrst observation, i.e. with 
the exchange rate on 03/13/2001. Figure 5.3 contains the graphs of the relative changes in 
exchange rates with respect to the US dollar for the same countries as in Figure 5.2[ We 
repeated our analysis for the relative changes (rescaled) in the exchange rates with respect to 
the US dollar, resulting in ^ 23,500 = 303 which corresponds to 05/24/2002. In the dehnition 
of ^ 23,500 we used M 2 — Mi ^ [T^/^/A 23 , 5 ooJ • The 90%, 95% and 99% conhdence intervals 
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Figure 5.4. The graphs of the log differences of the GDP/capita for 113 
countries between 1961 and 2012. 



are [292,330], [287,341] and [274,371]. Note that all these conhdence intervals contain 297 
which was obtained for the non-scaled exchange rates. Some of the graphs show a linear 
trend after the time of change instead of changing to an other constant mean. The limit 
distribution of the time of change is local, i.e. it is determined by the observations in the 
neighborhood of to- Replacing the linear trend with an average value close to to can justify 
the asymptotic validity of the conhdence intervals. Also, after the change point was found 
and the means of the corresponding segments were removed, the stationarity of the residuals 
could not be rejected. 

The exchange rates data and the scaled exchange rates contain further level shifts. Using 
the binary segmentation method one can divide the data into homogenous segments and 
provide conhdence intervals for the time of changes. Sato (2013) investigates the number 
and the location of the changes in daily log returns in 30 currency pairs between 04/01/2001 
and 30/12/2011 and points out that the study of the individual pairs might not hnd all the 
changes in the exchange rates mechanism. 


In the second example we compared the GDP/capita of = 113 countries. The data can be 
found at the website www://data.worldbank.org/indicator/NY.GDP.MKTP.CD. The data 
are recorded in current US dollar. We removed some countries from the data set due to large 
number of missing values, so we used = 113 panels with T = 51 covering the time period 
1961-2012. To achieve stationarity of the errors we transformed the data by taking log dif¬ 
ferences. The graphs of the log transformed GDP’s are exhibited in Figure 5.4[ We used the 
GUSUM test of Horvath and Huskova (2013) to test the stability of the means of the panels 
which was rejected at very high signihcance level. The estimator for the time of change is 
Ui 3 , 5 i = 19 which corresponds to 1979/1980. Applying the limit result in (4.6), [16,21] is 
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the asymptotic 90% and 95% confidence interval, while the 99% confidence interval is [12,22]. 


6. Conclusion 

We established the first and second order asymptotic properties of a CUSUM estimator 
of the time of change in the mean of panel data. Our results are derived under long run 
moment conditions, which serve to extend the asymptotic theory to a broader family of 
error processes than had been previously considered in the literature, and we provided an in 
depth study on how the rates of divergence of the sizes of changes and the common factor 
loadings are manifested in the asymptotic behavior of the test statistic. Our results were 
demonstrated with a Monte Carlo simulation study, and we considered application to two 
real data sets. 
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Appendix A. Proofs of Theorems I2.1H2.4I and Remarks 12.11 and 12.21 
Throughout the proofs c denotes unimportant constants whose values might change from 


line to line. Using (2.1) we have 


t 


where 


and 


Si{t) - —Si{T) = Qi{t) + 7iU (t) + 6ir{t), 


t T t T 

Qiit) = Vit) = ^Vs- 

S = 1 S=1 S = 1 S = 1 


r{t) = + > to}. 


Hence we have 


t 


(A.l) 

(A.2) 


S,{t) - -S,{T)j = Q2(t) + ^ ^ 2V{th,Q,{t) + 2r{t)6Mt) (A.3) 

-h 2V{t)r(t)'yi6i. 

Let 0<a<6'<l — a and dehne 

^ / t V 

tA,T(a) = argmaxLr„j<i<T-Lraj j . 

i=l ^ ' 


Lemma A.l. If Assumptions [AT - 2.5 are satisfied, then we have 


and 


lim PltjsfT = tisfT(a)} = 1 for all 0<a<9<l — a 

min(A,T)^oo ’ - ' 


^ 6 as min(T, N) ^ oo 


(A.4) 

(A.5) 
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Proof. It is easy to see that for every 0<a<6*<l — a 


We prove that 


— max r\t) ^ 6^(1-6?, 

T2 [Tai<t<T-[Tai ^ ^ ^ ' 

^ max r^(t) —)■ — 0)^, 

^ max r^(t) —)■ a^6^. 

T2 T-lTa\<t<T ^ ’ 


N 


sup 

i<t<r . , 
— — 2=1 


Y,QUt) ^ Op{NT). 


(A.6) 


Elementary arguments give 

EQ^t) < U,4t) + 2 + U,4T) 

and therefore by Assumption |2.3K i) we have 

N 

m^ax = 0{NT). 


2=1 


Let qi{u) = [QfiuT) — EQ‘^{uT))/T. Using Assumption |2.l[ i), for every 0<-u<n<lby 
the Rosenthal inequality (cf. Petrov (1995, p. 59)) we have with u = 

v/2 ' 


A 




i=l 


N 


N 


-QiiuW + '^E{qi{v) -qi{u)f 


2=1 


, 2=1 


By the Cauchy-Schwarz inequality we conclude for all 1 < s < f < T that 

E(QUt) - QMf < E{{QM - Q.(i!))^(|Qj(()| + |Qi(s)|)n 

< 4{E{Q,{t) - QM)y'^(EQf(t) + EQf(s)y'\ 
Using the dehnition of Qi(t) and Assumption 2.1[ii) we get that 


E(Q.(t)-Qi(s)y<2^{U,,4t-s) + 


t — s 
T 


u^AT) ^ • 


and similarly 
Also, 


EQt{t)<2^U,4t) + Ui,^{T)). 


{EQiit) - EQi{s)y < 2E{Qyt) - Qys)y{EQi{t) + EQiis)) 

2 


^ 8 I Ui^ 2 {t — s) + 


t — s 


U^,2{T)j (U,,2(t) + U,,2{s) + 2U,,2(T)) . 


Thus applying Assumption 2.3[ ii) we get that with some 0 < c < oo 

v/2 


N 


N 


X E{qi{v) — qi{u)y I < c|m — for all 0 < n < n < 1. 


2=1 













20 


LAJOS HORVATH, MARIE HUSKOVA, GREGORY RICE, AND JIA WANG 


Repeating the arguments used above we obtain that 

E\Qm - < E{mt) - + iq *( s )|)^} 

< 2^{E\Q,{t) - Q,{sTf/\E\Q,{tT + E\Q,{sTfl\ 


E\Q,{t) - Q,(s)|^" < 2^^^ - s) + 


t — s 
T 


2u 


u,.2Jt) ^, 


£| 0 .(«)R < 2"" (UisM + c/.,2„(r)) 

and 

\EQm - eqKs)\^^ < Emt) - Q.{sW{Emt)r+Ems)\^ 




t — s 
T 


U.,AT) (f/,,,(t) + f/,,,(s) + 2f/,,,(T)) 


resulting in 


N 


N 


E\qi{u) — qi{v)\'' < c\u — for all 0 < m, n < 1 


2 = 1 


with some c. Using Billingsley (1968, pp. 95 and 127) we conclude that the process 'Y^=\ Qii'Ej/N 
is tight in V[0, 1] and therefore (A.6) holds. 


The moment assumption in Assumption 2^ with the maximal inequality of Moritz et ah 
(1982) yields that i7(maxi<t<T |U(t)|)'' = 0(T'^) and therefore by Markov’s inequality we 
conclude 

max |U(t)| = Op(T^^^). 


i<t<r 


By (A.7) we get immediately that 


sup 


N 


l<t<T . , 
— — 2 = 1 


5^yr^(() = Op(i)Tr„,r. 


Following the proof of (A.6) we get 


sup 

l<t<T 


N 




and therefore by (A.7) 


sup 

i<t<r 


2=1 


N 


= OpmT^I^T'll. 




Similarly to (A.9) we have that 


sup 

i<t<r 


2=1 


N 


= Op(i)rr;;y 


T{i)5iQi{i) 


2 = 1 


op(i)T’-'M;,y. 


Using again Assumption |2.4| and the dehnition of r(f), one can easily verify that 


max 

i<t<r 


N 




2=1 


= Op(l)T3/2|S^,, 


(A.7) 


(A.8) 


(A.9) 


(A.IO) 


(A.ll) 
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It follows from (A.6)-(A.ll) that 


N 


T'^An,t 
1 


2=1 

N 


max5^(s.(()-yS.(T)) 4 


P 02, 


and 


Am 'T l<t<laT\ 

’ i=l 

1 




A ^. 2 , 


N 


T'^Amt T-\aT\<t<T ^ 

’ i=l 

which immediately implies Lemma [A. 1[ 




P. „.2n2 


□ 


According to (A.4), it is enough to consider the asymptotic behavior of 

iN,T = argmaxL„rj<t<[(i-a)rj {f^Jv(^) - f^A(^o)} , 
with any 0<a<6'<l — a<l, where UNit) is dehned in (4.2). It is easy to see that 
UN{t) — Uisfito) 


(A.12) 


N 


2=1 


+ Qlit) - Ql{to) + - V\to)) 

2'ji{Qi{t)V{t) - Qi{to)V{to)) + 26i{r{t)Qi{t) - r{to)Qi{to)) 

2'yi6i{V{t)r{t) - I/(to)r(to)) j- 


Lemma A.2. If Assumption \2.^ holds, then for all0<a<9<l — a there are 0 < ci, C 2 < 
oo such that 


-Ci\to - t\T < r^(f) - r^(fo) < -C2|to - t\T 


for all 1 < t < T. 


Proof. The result follows from Assumption 2.2 and the definition of r{t). 
Throughout Lemmas A.3 A.7 we assume that 1 < M < T. 


□ 


Lemma A.3. If Assumptions 2.1-A3 hold, then 


max 


|t-toT>M |t — to| 
Proof. We write 


N 




2 = 1 


= Op{N + + {NT/Mf/'^). (A.13) 


Qi W - QUIo) = 2Qi(to)(<5i(f) - Qi{io)) + {Qi{i) - Qi{to)f 


Using Assumption |2.3Ki) we get that 


N 


sup > 

l<t<T ““ 
— — 2 = 1 


W.it) - EQ^^{U)\ 
\t - tol 


= 0{N). 


(A.14) 
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Let Xi{t) = {Qi{t) - Qi{to)y - E{Qi{t) - Qi{to)y. Elementary arguments give 

1 

tn\ 

i=l 


P < max 


l<t<T \t -to\ 


N 




> X 


(A. 16) 


= P 


N 




i=l 


> x\t — to| for at least one 1 < t < T 


< 




l<t<T 


N 




2=1 


> x\t — to 


< {x\t-to\) 


l<t<T 


N 




2=1 


where in the last step we used Markov’s inequality. Let v = a/ 2 , where k is given in As¬ 
sumption 2.3[ ii). The processes Xi(t) are independent in so using the Rosenthal inequality 
(cf. Petrov (1995, p. 59)) we conclude with some c > 0, not depending on t. 


N 




i=l 


N 


N 


vl2 




2 = 1 


, 2=1 


2=1 


Assumptions 113 i) and |2.3| (ii) yield 

N 

'^E\xi{t)\'' < cN\t-to\^ 

and 

< cN-/'^\t - k\P 

Thus we conclude via ( A.15[ )-( p^.l7[ ) 


N 


vl2 


(A.16) 


(A.17) 


P < max 


l<i<T \t — to 


N 




2 = 1 


> X 


> < C > ^^ 


Choosing x = with a large enough c* we get that 

N 


max 


l<t<T |t - to| 


^XiY) 


2=1 




Let 




s=l 


It follows from the definition of Qi{t) that 


N 


- QM) - E[Qi{to){Qi{t) - Qi{to 


2=1 


(A. 18) 
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< 


N 


^{Qi{to){Si{to) - Si{t)) - E[Qi{to){Si{to) - Si{t))]) 


i=l 


\t - to 

T 


N 


5^(Q,(to)5,(T) - E[QMS,{T)]) 


2 = 1 


Using Assumption |2.lK i), and 2.3[ ii) with the Cauchy-Schwarz inequality we get that 

^ A \ N N 

J2Q‘MSi(T) = J]var(0.(i„)5i(r)) < J2iEQtMESf(T)y'^ = 0{NT^ 


var 


> *=i 
and therefore 


2 = 1 


2 = 1 


N 


j2mto)sm - E[QMs,{T)]) 


2=1 


= Op{nE^T). 


With (i(t) = Qi(to)ei^t — E[Qi(to)ei^t] we can write for 1 < t < to that 


N 


to N 




2 = 1 


S=t + 1 2=1 


By the Markov inequality we have 

1 


P < max 


l<t<to-M to — t 


to N 

S=t+1 2=1 


> xiNT/Mf/"^ 


< P i max max - 
logM<fc<log4o e*<£<e'“+l t 


to N 

S=tQ—£ 2=1 


< P< max - 
e'=<r<e'=+i £ 


log to 

^ E ^ 

fc=logM 


> x{NT/Mf/‘^ 
> x{NT/My/^ 

for at least one logM <k< log to 

> x{NT/Mf/‘^e^ 


to N 

s=to—£ *=1 


max 

e''<£<e*+i 


to N 

E E-i*!*) 

s=to—£ i=l 
log to 


< (a:(A^T/M)h2)-^ ^ ^-ku^ 

fc=log M 


max 

e'= <£<6*^+1 


to N 

E E 

S=tQ—i 2 = 1 


(A.19) 


Next we need a maximal inequality for double sum in the last term above. With Ci{s) = 
Ci(to — s + 1) we get that 


to N IN 

s=to—£+l i=l s=l *=1 
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By the independence of the processes Q^s) in i, Rosenthal’s inequality (cf. Petrov (1995, 
p.59)) implies that 


E 


N V 

1 = 1 s=u 


N 


< C 




2 = 1 




N / V 


v/2 ' 


We have via the Cauchy-Schwarz inequality 


^ Ec- 

\ S=U 

and therefore 


2=1 \s=U 


< {EQt{to)y'^E(J2ei 


(A.20) 


4 \ 1/2 


N / V 


N 


2=1 


^ {EQt(t„)y'^ (^(Ev. 

2=1 ^ ^S=U 

f N N / V 

E^«?(‘«)E^ Ev. 


4 \ 1/2 


1/2 


2 = 1 


2=1 \s=U 


Using that the Ci/s have mean zero and Assumption 2.3, we conclude 

N / V ^ 2 

<cNT\u-v\ 

2=1 \s=U 


Similarly, by the dehnition of Q we get 




<2^ I E 




+ 


E 


Qi{to) '^2 ' 


and by applications of the Cauchy-Schwarz inequality we have 




E 




< { EQl{t,)E 


< lE\Q,{to)rE 






1/2 


2u 


1/2 


resulting in 


N 


E^ 


2=1 


E<‘W 


< cNT'^>‘^\u-vV/‘^. 


Using the inequalities above, we get the upper bound for the moment in (A.20): 

V N 

E y/ y/ Ci{s) 


S=u 2=1 


< cN‘'/‘^T‘'/‘^\u - vr'‘l 


(A.21) 
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Applying the maximal inequality in Moritz et al. (1982) to (A.21) we conclude 

i N 


E max 

gfe<£<gfe+i 


S = 1 2=1 




Hence (|A.19|) implies that 
1 


P < max 


l<t<tQ-M tf) — t 


resulting in 


to N 

T. 

s=i+l i=l 


> x{NT/Mf/‘^ I < 


x^'^NT/NiyP 


11 J^V I2rp. 


•ujl 


k=log M 


C 


max - 

l<t<to-M Iq —t 

Similar arguments yield 

1 

max 


to+M<t<T to — t 

which completes the proof of the lemma. 


to N 

E Ec-i*) 

s=4+l i=l 

to N 

s=i + l 2=1 


= Op{l){NT/My/‘^. 


= Op{1){NT/MY/^, 


□ 


Lemma A.4. If Assumptions 2^- 2^ hold, then 

N 


max 


|t-to|>M \t - to\ 


'^Si{r{t)Qi{t) - r{to)Qi{to)) 


= Op{l){T^/^A](^^ + T(Ajv,r/M)i/2). 


2=1 

Proof. First we write 

5i{r{t)Qi{t) - r{to)Qi{to)) = r{t)6i{Qi{t) - QM) + SiQi{to){r{t) - r{to)). 
Applying the dehnition of r{t) with Assumptions |2.l| (i) and |2.3K i), we get 

N 


max 




\t-to\<T \t - to 


^ ^ diQiito 


2=1 




It follows from the dehnition of Qi{t) that for all i 


QM - QiY) = Zi{t,to) 


tn-t 


T 


if I <t <T, 


S=1 


where 


Zi{t,to) = < 


io 

ei,s, if 1 < ^ < to 

S=t + 1 

0, if t = to 

t 

- if to <t <T. 

S=to + l 
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Clearly, = 0{T on account of Assumptions 2.1'i) and 2.3'i) and 

therefore 


\r{t)\ 

|i-fo7<T |f - fo| 


max 


t-to 

T 


N 


2=1 S=1 




< 


N 


2=1 S=1 




Repeating the arguments used in (A.19), by Markov’s inequality we have 

N 

> x(A7v,T/Af)^^^ 

2 = 1 


P < max 


l<t<tQ-M tn — t 


y^^6iZi{t,to) 
1 


(A.22) 


< P < max max - 
logM<j<logT ei<£<ei+^ t 


N 


5jZj(fo—to 


2=1 




oo 

2 E ^ 


< 


j=log M 

(M/ 


max 

e^'<f<eJ+l 


N 


^6iZi{to - i, to 


2=1 


> xeYA^^T/My/^ 


X’' 


e max 


With Ci^s = we 


max 

eJ<r<e3+i 


j=logM 


N 


N 


YZi{to—i, to 


2=1 


J2YZYto-i,to) 


i=l 


= max 

eJ<£<eJ+l 


l+l N 

S = 1 2=1 


Using Assumptions 2.1[ i), [2^ ii) and 2.5[ ii) with Rosenthal’s inequality we conclude for all 
u > 2 that 


< 


V N 

E 

S=U 2=1 

since by the multinomial theorem 


N 


civ - uY'H Y,\YY + aYI U 2c(n - uYI’^aYI, (A.23) 


2 = 1 


TV 


5^ i«.r < 


2=1 


The maximal inequality of Moricz et ah (1982) and (A.23) imply that 


E max 

eJ<£<eJ+l 


N 


YYzyto-Yto) 


2=1 


and therefore by (A.23) we have 


P < max 


< c 


l<t<tQ —Af - t 

(M/ A7v,t 


N 


YYZiit,to 


2 = 1 

p/2 ^ 




> x(AAr^'r/M)^/^ 


(A.24) 


x^ 


j=log M 
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Thus we conclude that 


1 

max - 

l<t<to—M to — t 


N 

i=l 


and by similar arguments we have 


1 

max - 

tQ+M<t<T t — tg 


N 

i=l 


which also completes the proof of the lemma. 


Op{l){AN,T/My^‘^ 


Op{l){AN,T/My^‘^, 


□ 


Lemma A. 5. If Assumptions 2^- 2^ hold, then 

N 


1 


max 


\t-to\>M \t — to| 


J2invHt )- v ^ t ,)) 


2 = 1 


= Op(l)rjv,r (1 + (log(T/M))2/" + 


Proof We write \V\t) - V\to)\ < {V{t) - V{to)f + 2\V{to)\\V{t) - V{to)\. If 1 < t < tg, 
then 


\v(t) - V'lio)! < 


£ 

s=£+l 


+ 


\t - tg| 

T 


T 


s=l 


and therefore 


(r(<) - vm)Y < 4 ( ,Y + 4<t^ ( Y. ’'■) 

\s=t+l / \s=l / 


Thus we get from Assumption |2.4| that 
{V{t) - V{t,)f 


max 

\t-to\>M 


to -1 


= Op{l) + 2 max 


|t-to|>M {to — 


£ 

s= 4 +l 


Repeating the arguments used in (|A.19|) we get that 

io 


P < max 


to-t>M (to — tyA 

< p 


s=t+l 


logT 

A:=logM 
1 logT 


> X 

1 

logM<A:<logTe'=<tt<e*+i 
to 


max max 


to 

S = tQ—U-\-l 


> X 


max 


S=to—'ti+1 


> xe 


k/2 


< ^ V max 

X'^ e*<u<e'“+i 


fc=log M 


£ 

s=to—u-\-l 


Following the arguments used in the proofs of Lemmas |A.3| and |A.4| one can verify that 


E max 

e^<u<e^+^ 


io 

S = tQ — U-\-l 
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which implies that 


P < max 


to-t>M {to — 


to 




> X > < c 


log(T/M) 




Similar computations can be performed for t — to > M and thus we conclude 

to 


max 


s=t+l 


\tQ-t\>M {to - 

As in the proof of Lemma |A.3| we have that 
1 


= Op((iog(r/M))‘''*). 


sup 


to 


to—t>M to t 

The proof of the lemma is now complete. 


= Op{M and sup - 

t—to>M t to 


S=t0 + 1 


= Op{M-^/^). (A.25) 
□ 


Lemma A. 6. If Assumptions [^!T - 2.4 hold, then 

N 


max 




2 = 1 




\t-to\>M |t — to| 

Proof. We write 

Q^{t)V{t) - Qi{to)V{to) = V{t){Q,{t) - Qi{to)) + {V{t) - V{to))Q^{to). 


Assumption 2.4 implies that 


max |l/(t)| = Op(Ti/ 2 N 

l<t<T 

and by the arguments used in the proof of Lemma |A.4| one can show that 

1 


max 


\t-to\>M \t - fo| 

Similar arguments yield 


N 


- QM) 


2=1 


= Op{l)M-y^T]l^. 


max 


\V(() - 1^(4)! 


\t-to\>M \t — to 


N 




2=1 


Op{l)T^/‘^M-^yT^J^%. 


□ 


Lemma A. 7. If A ssumptions \Kl - 2.4 hold, then 

N 


max 


\t-to\>M \t - to| 


- v{to)r{to)) 


2=1 


= Op(l)(T ^/2 ^ TM-y^) ISat^pI 


Proof. Since V{t)r{f) — V{to)r{to) = V{to){r{t) — r{to)) +r{t){V{t) — V{to)), Lemma [AT 
follows from 

v(t„)PPPP> =Op(i)r‘/^ 


max 

\t—to\>M 


and 


max 

\t-to\>M 


r{t) 


t - to 

V{t) - V{to) 


t-to 


= Op{l)TM-y^. 


(A.26) 

(A.27) 
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The claim in (A.26) is an immediate consequence of the dehnition of r{t) and Assumption 
^ while (I A. 271) is proven in (|A.25|). □ 


Proof of Theorem 2.1 Under assumptions (2.3) and (2.4) we use Lemmas A.2-A.7 with 

M = 1. □ 


Proof of Remark 2.1. The proof of this remark follows Bai (2010) closely. We use (A.3). 
Since T is hxed, 

N 


i=l 


and by Assumption 2.3[ i) and Markov’s inequality we have 

N 

= Op{N). 

~ ~ i = l 

By the Cauchy-Schwarz inequality and Assumption 2.3[i) we conclude 


N 




and therefore 


Similar arguments give 


and 


1=1 

N 




max 

l<t<T 




2 = 1 


= 0p{i)Ty^. 


max 

l<t<T 


max 

i<t<r 


N 




2=1 


N 


it)y{t)'^Sai 


= OpiAy^) 


= Op(|E|). 


2 = 1 


The hnal term coming from (A.3) to consider is T,fLir‘^(t)6f = Ar^(t). Under the conditions 
of the remark, this is the asymptotically dominating term which has a unique maximum at 
to- Hence Remark 2^ is proven. □ 

Proof of Remark \2.^ Let 


= 


(t(T-t))i/2 




where Qiif) and U(t) are dehned in (A.l). We note that due to the assumption that the et^s 
and rjt are sequences of uncorrelated random variables we get that 

= + (A.28) 

and 

var(/2(t))<C,(Ueto + 7f) (A.29) 

with some constant Ci. We write 

N / , N 2 1 


2=1 




T 


t{T - T) 


— 'kLl,N{i) + '^ 2 , a (^) + 'kLz,N{l)i 
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with 


n 2 ,Nit) = fKt) and H 3 ,n(t) = 


N 


N 


6ir{t) 


t{T - T) 


2=1 


2=1 


(t(T-t))V 2 ’ 


where r{t),t = 1, 2,..., T is dehned in (A.2). We show that for all t ^ to 


lim P{'Hl,Ar(^o) — 'Hl,N{t) < 'H 2 ,N{t) — 'H 2 ,N{to) + 'H 2 ,N{t) — 'H 2 ,N{to)} — 0, (A.30) 


N^OO 


which immediately implies Remark |2.2[ We note that with some 6*2 > 0 we have that 
H i,a(^o) ^ C 2 A for all t ^ to. By the independence of the processes Qi{t')^ 1 ^ 

i < N and V(t) we conclude 


N 


E{n2,N{t)-'H2,N{to)f = Y^ 


2 = 1 


:QKt) 


N 


4Y,E 


t{T-tr^^ ' to{T-to) 
1 


2 = 1 


f(T - 1) 


Qi{t)%v{t) - 


QKto] 

1 


~ to) 


QM'yiVit) 


E 


V\t) V\to) 


XT - t) to{T - to) 
= 0{N + T + X) 


and therefore 


max^ \n 2 At) - %,iv(to))| = Op{l) + p) 


Similarly, 

m^ax |'R3,7v(f) - 'R3,Ar(fo))| = Op{l) + |S|) = Op(l) (A^/2 + , 

since |S| < completing the proof of ( |A.30 ). 


□ 


Lemma A. 8 . VBe assume that Assumptions \2. l]\2X hold, and |s| < 00 . Then, as N,T -^00 
we have that 


AN,T\iN,T — ^o| — Op{l). 

Proof. By Lemma A.l it is enough to prove that for all 0 < a < 0 

A\i]y,Tioi) — to| = Op{l) 


(A.31) 


(A.32) 


Under Assumption |2.6[ i) we choose M = C/Apf^p, where (7 > 0 is a constant. Using Lemmas 
A .2 |A.3| |A.5[jA.7] and Assumption |2.5[ ii) we obtain that 


N 


A«,T(r^(i) - r=(i„)) + - QlM) + r«.T(v"=(i) - v^M) 


(A.33) 


i=l 


N 


N 


‘^YT(Qi(t)V{t) - Qi{to)V{to)) + 2^7ihi(U(f)r(f) -V{to)r{to)) 


2=1 
.2 


2=1 


= Apf^T^it) — m{to)){^ + Op{l)) uniformly on \t — to\>M 
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for all M. Also, by Lemmas A.2 and A.4 we obtain that 


N 


lim lim inf P 

C^oc N,T^oo |t-to|>C/Ajv,T tiT 


sup ^ - r'^ito)) + ‘26i{r{t)Qi{t) - r{to)Qi{to))] < 0 


Hence Lemma A .8 is established nnder Assnmption |2.6[i) and | 5 | < oo 


1 . 

□ 


Lemma A.9. 

we have that 


We assume that Assumptions 2.1-2.6 hold, and |s| < oo. Then, as N,T —)■ oo 


snp 

\t-to\<C/A 


f - r^{t)) - 20(1 - e)Age{t - to) 

i=l 


o(l), 


(A.34) 


snp 

\t-tQ\<C/A 


^ ^ di{Qi{t)r{t) - Qi{to)r{to)) 


— op(l), 


N 

+ e{l-e)J2^^is^it)-s^{to)) 

i=l 


(A.35) 


snp 

\t-to\<C/A 


^ ^ lA{V (t)r(t) - V(to)r(to)) + 0(1 - 6)En,t(V( t) - V(to)) 


snp 

|t-to|<C/A 


N 


T 


-«?(*»)) 


snp 

lt-tol<C/A 


i=l 

N 


T 




— Op(l), 
= Op(l) 


2 = 1 


and 


snp 

\t-to\<C/A 


N 


^ii{Qi{t)y{t) - Qi{to)V(to)) 


2 = 1 


- 0p(l), 


for all C > 0, where A = Aat^t and Si(-) is defined in (A.18). 
Proof. First we note 


1 

T 


- r^^ifo)) = - r(to)) 


2 = 1 


2 = 1 


N 


T 


'^d‘l(r(t) - r(to)f 


2=1 


Using the definition of r(t) and Assnmption 2.5[i) we conclnde 


snp 

|t-to|<C/A 


N 


T 


'^d‘i(r(t) - r(to)f 


2 = 1 


= 0(1/(TA)) = 0(1) 


op(l), 

(A.36) 

(A.37) 

(A.38) 


(A.39) 
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and 


sup 

|i-tol<G/A 


- r{to)) - 20(1 - e)Agg{t - to) 

i=l 


= o(l), 


completing the proof of ( A.34[ ). 

Similarly, 

N 

Si{Qi{t)r{t) - Qi{to)r{to)) = r(to) 5i(Qi(t) - Qi(to)) + 5iQi(t)(r(t) - r(to)) 


N 


N 


1=1 


2=1 


2=1 


and 


N 


N 


-Q*(to)) = 5.(5,(t) -5,(to)) + ^^5^5.5.(T). 


N 


2=1 


2=1 


T 


2 = 1 


Computing the variance of 5.5, (T) we get 


sup 

\t-to\<C/A 


N 




2=1 


= Op(i/(rA)>/2) = op(i) 


by Assumption |2.5K i), so (A.35) is proven. 

Clearly, 

V (t)r(t) - V (to)r(to) = (V (t) - V (to))r(t) + V (to)(r(t) - r(to)). 


By Assumption 2.4 we get that C(to) = Op(T^^^) and therefore 


sup |C(to)(r(t) -r(to))| = Op(T^''^/A). 

lt-tol<C/A 


We note that for all to < t < to + C/A 

|rw-r(((,)|< 


E ■■ 

s=to+l 


+ 


l^-^o 

T 


E 

s=i 


Vs 


and by Assumption 2.4 we have that | Yll=to+i Vs\ = Cp(l/A^/^) and the process A^/^ Vsy 0 < 

M < 1 is tight in P[0, C]. Thus by stationarity we get 


sup 

to<t<C/A 


ris =Op(l/A^/2) 

and similar arguments can be used on to — C/A < t < to. We conclude that 


sup |C(t) - f/(to)I = Op(l/AV2 + T-1/VA) = Op(l/Ai/2) 
\t-to\<C/A 


(A.40) 


completing the proof of (A.36) on account |s| < oo. 

With '0,(t) = Q‘^{t) - Qjito) - E{Q‘l{t) - Qf{to)) we can write 

N 


sup 

\t-to\<C/A 


T 


EW?(«) - QlM) 

1 


2=1 


< sup — 


\t-to\<C/A 


T 


N 


J2(EQKt) - EQUto)) 


2 = 1 


+ sup — 


|t-to|<G/A 


T 


N 




2=1 
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We obtain from the proof of Lemma A.l that 


1 

sup — 

\t-to\<C/A T 


N 


X^(EQUt) - EQfM) 



For every t G [to — ^F/A, to + C'/A] we have that 


N 


N 




i=l 


2^2 


i=l 


and 


0 ( 1 ). 


sup E'ljjfit) < 9 sup {Q- (t)((5i(t) - Qi{to)f + Q- (to)(Qi(t) - Qi{to))^} 

\t-to\<C/A \t-to\<C/A 

= 0(1) I sup (EQ^(t))^/2 sup (E(gi(t) - Qi(to))^)^/^ , 

[|t-to|<C/A \t-to\<C/A J 

= 0{T/A). 


Next we show that a/A/ (NT) J2iLi A) is tight in V[—C,C]. Using Rosenthal’s in¬ 
equality (cf. Petrov (1995, p. 59)) we obtain that 


E 


N 




i=l 


k/2 


< 


N 

cl '^E\'ijji{t) -V'i(s)r^^+ I '^Ei'ipiit) -V’i(s))' 


N 


k/4 ' 


2 = 1 


. 2 = 1 


with some constant c. It is easy to see that 

- '0i(s)| < {|gi(t)(gi(t) - gi(s))| + |gi(s)(gi(t) - gi(s))| + i^g-(t) - (s)|} 


and 

|ug2 (t)-ug2(5)| <c|t-s|. 

By the Cauchy-Schwarz inequality and Assumption|2.3|[ii) we have for all t, s G [to—O/A, to+ 
C/A] 

E{Q,{t){Qi{t) - Q,{s)) f < {EQt{t)E{Q,{t) - Qi{s))y/^ < cT\t - s\ 


and therefore 

E{iJi{t) - i/’i(s))^ < cT|t - s| 

where c is a constant. Also, for k of Assumption 2.3[ ii) we have 

Emt){Q,{t) - g.(s))r/' < {U|g,(t)ru|g,(t) - q^{sTY/^ < c{u,,YmA\t - 

with some constant c. Thus we get via Assumption |2.3[ii) that 


A 

Wf 


1/2 N 


'Y^{i)i{u/A) -'ipiY/A)) 


2=1 


k/2 ' 


< c\u 


|k/4 


establishing tightness by Billingsley (1968, pp. 95 and 127). This also completes the proof 
of (lAATl). 


Following the arguments in the proof of (A.36) one can show that 


sup |U^(t) - U^(to)| < 2 sup |U(t)-U(to)| sup \V{t)\ = Op{T^^^/A^^Y^ 

\t-to\<C/A \t-tQ\<C/A \t-to\<C/A 
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and therefore ( A.38[ ) follows from Assumption 2.5[ ii). To prove (A.39) we first write 

Q^{t)v{t) - QMV{to) = V{t){Q,{t) - QM) + Q^{to){V{t) - V{to)). 


Repeating the arguments used in the proof of (A.37) we obtain that 

N 


T 


sup 

|t-to|<C/A 




2 = 1 


= Op(l) 


r \ 

i N,T 

~fA 


= Op(l) 


via applying Assumption 2.5 ii) and TA —)• oo. □ 

Let 

N 

Rn,t{u) = ''^^5iSi{u/ u>0, 

2 = 1 

where iSj(-) is defined in ( |A.18[ ). 

Lemma A.10. If Assumptions \Kl\ \2.^ ( |2.8| ) and ( |2.10[ ) hold, then we have 

Rn,tW —t aW{u), 

for all C > 0, where W{u) stands for a Wiener process. 

Proof. For the sake of notational simplicity we write A = A^r^p. Let 0 = Mq < < "^2 < 

... < Uk < C and ai, ..., a^. Under (2.10) we write 

k N k 

<^iiRN,T{ui) — Rn,t{ui-i)) = A) — Si{ui_i/A)). 


e=i 


i=i e=i 


Using Assumptions 2.1, 2.3, 2.7(i) and (2.10), we get that 


N 


^ E '^ai6i{Si{ue/A) - Si{ue-i/A)) = a'^'^aj{ue - ue-i){l + o(l)). 


i=l 


e=i 


£=1 


Also, Assumptions 2.1 (ii) and 2.7 li) imply 


N 




2=1 


Si ^ ai{Si{ui/A) - Si{ui_i/A)) 


1=1 
N 

2=1 


(A.41) 


max Ui f{\uc - Wf_i|/A) 

i<e<k 


N 


<c 


2=1 


So using Lyapunov’s theorem (cf. Petrov (1995, p. 154)) we conclude via (2.8) that 

k k 

C£t{RN,T{ut) — Rn,t{u£-i)) —)■ a ^^a£(lU(M£) — W{uii-i)), 

£=1 1=1 

where W stands for a Wiener process. Applying the Cramer-Wold theorem (cf. Billingsley 
(1968, p. 49)) we obtain that the hnite dimensional distributions of Rn,t{u) converge to 
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that of aW(u). Next we show that Rjsi^Tiu) is tight in P[0,C]. Following the arguments in 
( A.41[), Rosenthal’s inequality yields for all 0 < < C 


E\Rj^^'j'(u) — Rmt(v 


N,Ty 


N 


N 


f/2' 


< cA -"/2 J ^ \S,fE\Si{u/A)-Si{v/A)f + Y,6lE{S,{ul A) - S,{v/A)f 


i=l 


, i=l 


N 


N 


r/2' 


< d E- I-I/A) + E'*?L.2(|f - !>I/a) 


2=1 


N 


<c<^ A-"/2^i5,r + i 


\u — v\ 


. 2=1 


f/2 


2=1 


< c\u — 


on account of Assumption 2.7 ii) and (2.8). The tightness now follows from Billingsley (1968, 
p. 127). □ 


Lemma A.11. If Assumptions 2.1, 2.3, 2.1 and (2.8) hold, then for all integers 0 < ti < 
t 2 < ... < tx we have that 

N 




V 




„ 2=1 


where the Gaussian process &{t),t = 0, ±1, ±2 .... is defined in Theorem 2.2 


Proof. We repeat the hrst half of the proof of Lemma A.10 The result of Lemma A.11 
follows from (A.41) and Lyapunov’s central limit theorem due to assumption (2.8). 


□ 


Proof of Theorem 2.2. Let A = A^v.t- R follows from Assumption |2.1[ ii) and Lemma A.10 
that for all C* > 0 


A 


Y,d,{S{RRulA)-S{R)) 


v[-cfr\ 


aW{u), 


(A.42) 


2=1 


where hF(-u),— cxd < m < cxd is a two sided Wiener process. Also, Ago{u/A) = ge{u) and 

(A.43) 


since 5 = 0 by (A.40) we have that 


sup |R(f) - l/(to)| = op(l)- 

|t-to|<A 

By Lemma [A.9| we conclude that for all C > 0 

^ (RA(to + u/A) - U^ito)) 20(1 - e){aW{u) - ge{u)). 


(A.44) 


By the continuous mapping theorem we conclude from (A.44) that for all C 


(D 

argmaX|i_i^l<c-/A (^^(^0 + u/A) - f/7v(fo)) —t argmaX|„|<(^(crlF(M) - ge{u)). (A.45) 

According to the law of iterated logarithm, we have that 

Jim argmaX|^l<(^(crhF(M) — 5fe(M)) —t argmax„((ThF(M) — 5 fe(M)) 


a.s. 


(A.46) 
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Now (2.11) follows from Lemma A.8, (A.45) and (A.46). 

It follows from Assumption 2.1[ii), (A.43) and Lemmas A.9, A.11 that for every integer 
C>0 


y {UN{to + t) - UN{to)), t = 0, ±1, ±2 ,..., ±C 


(A.47) 


V 


{26(1 - 0)(0(t) - dgeit)), t = 0, ±1, ±2 ,..., iC} . 

Observing that u{t,t) = 0{t), the normality of 0(f) with the Borel-Cantelli lemma yields 
that lim|i|_j.oo = 0, and therefore 


Jim argmax|^|<(^(0(f) — dgeit)) = argmax^(0(f) — dgeit)) a.s. 


The proof of (2.13) is now completed via Lemma A.S 


□ 


tion 


Proof of The orem\2.^ It follows from Assumption 2.8 that the Wiener processes in Assump- 
2.9 and (|A.42) are independent. Hence Lemma A.9 yields for all C* > 0 that 


— (t/Ar(fo + m/A) — f/Ar(fo)) 




2e{l-e){a‘^+s^f/^W{u)-ge{u)), 


where W{u)^ —oo < m < cxd is a two-sided Wiener process. Arguments used in ( A.45| ) and 
(A.46) could be repeated to hnish the proof of (2.15). 


Referring again to Assumption 2.8 it is immediate that the Gaussian process 0(f) and V(f) 


are independent. So applying Lemma A.9 we replace (A.47) with 


^ {UN{to + t) - GAr(fo)) , f = 0, ±1, ±2, . . . , ±C 


V 


{20(1 - 0)(0(f) + - ^geit)), t = 0, ±1, ±2 ,..., ±C} 


for any G > 0. Observing that V(f)/f —)■ 0 a.s. we need only minor modihcations of the 
proof of (2.13) to complete the proof of (2.16). □ 


Lemma A.12. We assume that Assumptions 2.1-2.6 hold, and |s| = oo. Then, as N,T —>■ 
oo we have that 

|f7v,r — fo| = Op{M]y^T), (A.48) 

where = (Sjv,T/AAr^'r)^. 


Proof. By Lemma [A. 1 it is enough to prove that for all 0 < a < 0 

|^v,t(d) — ^oI = Op{Mn^t)- 


The result follows from Lemmas A.2-A.7 with M = (S 7 v,r/AAr,r)^- 


(A.49) 

□ 


Proof of Theorem 2.f. Let M = ■ Since lS.]si;r is bounded, by (2.17) we have 


A 


—)■ oo. 


N,T 


Following the proof of Lemma A.9 one can show that for all G > 0 

N 


sup 

|i-tol<GM 


di{Qi{t)r{f) - Qi{to)r{to)) 


2=1 


= C'p(l), 


(A.50) 


(A.51) 
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sup 

\t-tQ\<CM 


- QlM) 


sup 

\t-to\<CM 


i=l 

N 


T 




= Op(l), 
= Op(l) 


2=1 


and 


sup 

\t-tQ\<CM 


N 


^'yi{Qi{t)v{t) - Qi{to)v{to)) 


2=1 


= op(l)- 


(A.52) 

(A.53) 

(A.54) 


It follows from Assumptions |2.4| and |2.9| that for all C > 0 


— (AAr,r(r^(fo + uM) - r^(to)) + 2En,t{V{ to + uM)r{to + uM) - V(to)r(to))} 


T 


V[-C,C] 


2d(l-d)(W(u)-gs(u)), 


where W{u)^ —oo < u < oo denotes a two-sided Wiener process. Using now (A.51)-(A.54) 


^{UN{to + uM)-UN{to)) 29{l-9){W{u)-g0{u)). 


Arguments used in (A.45) and (|A.46) could be repeated to complete the proof of Theorem 

El □ 


Appendix B. Proof of Theorem 14.11 

For the sake of brevity we use t, A and S for tN,T, Ajsi^t and respectively . We start 

with the proof of (4.4). Let M = M{N,T) be a sequence satisfying 

M —)■ oo, M/T —)■ 0, (B.l) 

if the conditions of Theorem 12.21 or 12.31 hold and 

M ^ oo, M/ min(T, T?/A^) 0 


(B.2) 


under the assumptions of Theorem 2.4 First we show that for all M satisfying M/T —)■ 0 
we have that 


sup 

|n|<M 


N 

E 


E •’^7 


A \to + u 

1=1 \ ^ l<t<to+-u 


T - {to + u) 


E 


tQ-\-U<t<T 


— op(l). (B.3) 


Using (|2.1[) we have for all 0 < n < M 


1 


to + u 


E — 

l<t<tQ-\-U 

1 


^ ~ (^0 + n) 

E “ 


E -T,. 

tQ+u<t<T 

li 


toAu 

E 




to + u 


Y. vt 




E 


to + u T — {to + u) ^ T — {to + u) 

^ l<t<to+u ^ ^ ’ T-{tQ^u)<t<T ^ ^ ’ r-po+n)<4<r 




Since M/T —)■ 0, we have 


sup 

0<2i<M 


^0 


to + u 


0 . 


(B.4) 
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Applying Assumption 2.4 
we obtained for all z > 0 




P 


r 


I z 

T A o<u<M 


{ > ^ > = p I sup { V 

\l<t<tQ+u J J |^0<n<M 


I,,] >(;rM/r)*-'4 (B.5) 


< 


\T Aj „<„<M y<s^+„ J 


account of Assumption 
fo oTi^ lo Ql -^g conclude 


/ F \ 

.0(1) -0 


on account ol 
^i) and|^ 


2.5 


' ' \TAJ 

Following the proof of (B.5) but now using Assumptions 


"-1 \l<t<to 


•{^i:(e4>4<z^i 

by Assumption 2.5(i). The stationarity in Assumption 2.1[ ii) yields 

/ ,r / V 9 ^ 


(B.6) 


P 


1 


- i=l \to+l<t<to+u 


ei,t > z 


(B.7) 


P 


Vi£ 


ei,t > z 


M I 1 ^ 
u=l i=l 


^ NM 
= O 1 Hatw ^ 0 
^ ’ zT^ A 

by Assumption |2.5K i) and the assumption that M/T —)■ 0. Putting together (B.6) and (B.7) 
we obtain that 

X s 9 


E E O:,. 

i=l \l<t<u 


> Z 




Following the proofs of (B.5) and (|B.8|) one can prove that 


1 Y / 

0<U<M A ^ (^T - (to + u) 


and 


Vt 


N 


oSl<M A £ (^T - (to + u) 




op(l) 


op(l). 


(B.8) 


(B.9) 


(B.IO) 
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The result in (B.3) now follows from (B.4), (B.5) and (B.8)-(B.10). Since under our condi¬ 
tions 

\i-to\/M = op{l), (B.ll) 


the proof of (4.4) is complete. 


By (A.3) we have 


UN{t + n) - UnH) = U'^^\v) -1- ... + 


Ti^)l 


where 


N 


N 


i=l i=l 

N N 

uS\v) = + t.) - V\i)), ui^\v) = 2^7i((3i(f + v)V{i + v)- QSmi)). 


i=l 

N 


i=l 


Ui ^\ v ) = 2^6 i { r{i + v ) Qi{i + v ) - r { i ) Qi { i )), 

i=l 

and 

N 

Ui^\v) = 2'^Si'yi{r{i + v)V{i + v) -r{i)V{i)). 


2 = 1 


It follows from the proofs of Lemmas | A.3 , A.5, A.6 and A.9 that for all M satisfying (B.l) 
we have 

1 r ... .2 


L ,|^2 =Opil 

\v\<M N,T ^ ^ 


^N,T), 


where fpi^p is defined in (4.3). Applying now Lemmas A.4, A.7 and A.9 we conclude 


sup 


and 


where 


where 


I 

|<M I 


sup 


v\r 


~ — Op{'^N,t) 


N,T 


^ - uPiv) ' = Op(S„,3.), 




N 


N 


Uj^\v) = 2r{to) ^ SiAi{v; i) and UP{v) = 2r{to) ^ 6aiB{v; i), 


2=1 


2=1 


/' t-\-v 

E 

s=t+l 


ei,s, V > 0 


Ai{v;t) = < 


0 , 


n = 0 


t-i 

E' 

V. s=t-\-v 


V < 0, 
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and 


f t-\-v 


S=t+1 

B{v; t) = I 0, 


^ Vs, v>0 
V = 0 


t-i 


^ Vs, V < 0, 


s=t-\-v 


Using (B.l) and (B.ll) one can verify that 

1 


and 


where 


snp 

\v\<M 


snp 


\v\fN,T 

1 


u'n\v) - £/y(to 


( 9 ), 


— Op[^N,T) 


\v\<M \v\fN,T 


=Op(S»,t). 


r(l0) 


N 


N 


= 2r(to)'^SiAi{v,to) and = 2r{to)'^6aiB{v,to). 


i=l 


2 = 1 


Let m = m{N,T) < M and m —)■ oo, as niin(A^, T) —>■ oo. Using Assnmptions 2.3 and 2.4 
we get 


snp 

m<\v\<M 


4|n|r2(to) 




( 10 ), 


' N 


N 


^N,T 


and 


snp 


|i,|<m4|n|r2(to)S7v,T 


We claim that 


snp 

m<\v\<M 


4:\v\A{to)EN^T 


{uS\v) + uii’‘\v))' 


— o{En,t)- 


Op{l). 

(B.12) 

1 = Op(l). 

(B,13) 


The statement in (B.13) is a nniform weak law of large nnmbers, so we can repeat the proof 


of (|B.4|) to prove it. Namely, dne to stationarity, it follows from Assnmptions 2.3 and 2.4 

,, M] that 

(U’U) + U‘°A’)) 


that for every v G [—M ,..., —m, m,..., M] that 

1 


A\v\A{to)EN^T 


4 1 , 


(B.14) 


Now (B.13) follows from ( B.14[ ) if (u^\v) + /{\v\A{to)Eisf^T),m < |n| < M is 

tight. The tightness can be proven along the lines of the proofs of Lemmas A.4 and A.7 
Onr argnments show that 


snp 


4|y|f2 


N,T 


{UN{i + v) - UN{i) - AN,T{r^{i + v) - r^(t)))' 


(B.15) 


\v\r{to 


uS\v) + Uriv) 


( 10 ), 


— Op[^N,T> 


Pntting together (B.12), (B.13) and (B.15), the resnlt in (4.5) follows. 








































ESTIMATION OF THE TIME OF CHANGE IN PANEL DATA 


41 


Lajos Horvath, Department of Mathematics, University of Utah, Salt Lake City, UT 84112- 

0090 USA 

Marie Huskova, Department of Probability and Mathematical Statistics, Charles Univer¬ 
sity, SOKOLOVSKA 83, CZ-18600 Praha, Czech Republic 

Gregory Rice, Department of Mathematics, University of Utah, Salt Lake City, UT 84112- 

0090 USA 

JiA Wang, Department of Mathematics, University of Utah, Salt Lake City, UT 84112-0090 
USA 



